Magnetic phase diagram of a frustrated ferrimagnetic ladder: Relation to the 

one-dimensional boson Hubbard model 
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We study the magnetic phase diagram of two coupled mixed-spin (1, |) Heisenberg chains as 
a function of the frustration parameter related to diagonal exchange couplings. The analysis is 
performed by using spin- wave series and exact numerical diagonalization techniques. The obtained 
phase diagram - containing the Luttinger liquid phase, the plateau phase with a magnetization per 
rung M = 1/2, and the fully polarized phase - is closely related to the generic {J/U, IJ./U) phase 
diagram of the one-dimensional boson Hubbard model. 
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Some special classes of quantum spin ladders composed 
of different types of spins and/or regular mixtures of 
ferromagnetic and antiferromagnetic exchange bonds ex- 
hibit a number of novel quantum spin phases and unusual 
thermodynamic properties;:^'^ Moreover, such systems 
are appropriate for studying magnetic quantum critical 
points in one spatial dimension (ID), since many of them 
sustain magnetically ordered ground states. Possible re- 
alizations of such ladder structures can be based, for in- 
stance, on the synthesized quasi-lD bimetallic molecular 
magnets 1^ 

In a recent paper we considered the role of geometric 
frustration in a system of two coupled mixed-spin (1,5) 
Heisenberg chainsj^ It has been argued that a relatively 
moderate strength of the magnetic frustration produced 
by the nearest-neighbor diagonal bonds was able to de- 
stroy the ferrimagnetic phase and to stabilize the Lut- 
tinger liquid phase. The quantum phase transition be- 
tween these two phases is smooth and is realized through 
an intermediate canted spin phase. This latter phase is 
characterized by a net ferromagnetic moment per rung, 
< M < 0.5, and power-law transverse spin-spin corre- 
lations. In this Brief Report we show that the magnetic 
phase diagram of the discussed model as a function of the 
frustration parameter is closely related to the {J/U, n/U) 
phase diagram of the boson Hubbard model in periodic 
ID potentials,^ where /^, J, and U are, respectively, the 
chemical potential, the hopping strength constant, and 
the on-site repulsion potential. 

The relevant spin Hamiltonian reads 



W — 2^ [ Jl (Si„ • S2„+i + S2„ • Si„+i) -I- J_lSi„ • S2„ 
ri=l 
+ J2 (Sl„ • Si„+i + S2„ • S2„+i) - gfiBH{si^^ + S2^J] , 

(1) 

where L, g, fiB, and H are, respectively, the number of 
rungs, the gyromagnetic ratio, the Bohr magneton, and 
the external magnetic field applied in the z direction. 
The spin operators Si„ and S2„ (characterized by the 



quantum spin numbers Si > S2 and the rung index n) 
are alternatively distributed on the ladder sites. The 
parameters Ji , Jj_ , and J2 are positive coupling constants 
of the nearest-neighbor horizontal, vertical, and diagonal 
exchange bonds, respectively (see Fig. 1 in Ref. |5l). We 
set the energy and length scales by Ji = 1 and ao = 1, 
where ao is the lattice spacing. For simplicity, in the 
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FIG. 1: Phase diagram of the model Q in the (J2, h) plane. 
The phase boundary hs = hs{J2) of the fully polarized phase 
[M — 3/2) is exact. The dashed curve abed marks the area 
occupied by the M = 1/2 plateau phase, as obtained in a 
linear spin- wave approximation. The boundary position of 
this phase, as obtained from exact numerical diagonaliza- 
tion of periodic clusters, is traced by squares (L = 10) and 
crosses (L = 12). The dotted curve shows the phase bound- 
ary where Aq (/i) — [see Eq. Q], as obtained from the 
first-order spin-wave theory. The dashed curve starting from 
the point h and going upwards left as well as the piece fee 
represents in a linear spin-wave approximation the boundary 
where Ai" (/i) = 0. The canted spin phase exists for /i = in 
the interval 0.342 < J2 < Jc, where Jc = 0.399 is the phase 
transition point from the canted phase to the Luttinger liquid 
phase. The rest of the area is occupied by the Luttinger liquid 
phase. 



following we consider the case J± = 1 and set h — gfisH. 

The magnetic phase diagram of the model (Q) in the 
( J2, h) plane is presented in Fig. Q] The diagram contains 
the Luttinger liquid phase, the plateau magnetic phase 
with a quantized magnetic moment per rung M — 1/2 
and the fully polarized phase with M = 3/2. In addition, 
there is a canted magnetic phase which exists only at zero 
magnetic field (in the interval 0.342 < J2 < Jc = 0.399) 
and for /i > merges into the Luttinger liquid phase. 
Some properties of the latter phase at h = have al- 
ready been studied in Ref. |^ by using a conformal-field- 
theory analysis of the numerical exact-diagonalization 
(ED) data for the ground-state energy and the energies of 
some low-lying excited states in finite periodic systems. 
As is well known, the external magnetic field varies the 
parameters of the Luttinger liquid phase without chang- 
ing its basic properties up to a saturation field (hg), where 
the system becomes fully polarized,^ 

In this connection, let us first consider the critical field 
hs defining the phase boundary of the fully polarized 
M = 3/2 phase. The boundary hs = hs{J2) can be ob- 
tained by examining the instabilities of the one-magnon 
excitations above the ferromagnetic vacuum. We assume, 
as generally accepted, that the multimagnon excitations 
become unstable at a weaker magnetic field as compared 
to the one-magnon excitations. As a matter of fact, the 
latter assumption is confirmed by our ED data, which 
exactly reproduce the analytical result for hs presented 
below. A straightforward calculation results in the fol- 
lowing exact expression for the one-magnon dispersion 
relations in the fully polarized phase 
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where the wave vector k runs in the lattice Brillouin zone 
— TT < A; < TT. Equation Q describes two folding ferro- 
magnetic (FM) branches of one-magnon excitations. The 
phase boundary hs = /is( J2) is defined by the expressions 



hs = ^ (0 < J2 < 0.5) 
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where the first and second lines correspond, respectively, 
to instabilities of w^^^ (fc) at the wave vectors k — and 
TT. For h < hs{J2), the modes around k = Q and -k 
become unstable and destroy the fully polarized state. 
As discussed in the following text, the point E in Fig. 
n is a crossover point at which both minima in cij^^'(fc) 
reach the ground state. 

Let us now address the M = 1/2 plateau phase. The 
presence of quantum fluctuations in the ferrimagnetic 



state does not allow exact considerations, so that we shall 
rely on a qualitative analysis based on the spin- wave the- 
ory (SWT) supplemented by numerical ED data for finite 
periodic systems containing up to L = 12 rungs. In a lin- 
ear SWT approximation, the magnon excitations above 
the classical ferrimagnetic ground state are described by 
the following dispersion relations: 
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In the last equation, w'+-'(fc) is a FM branch, whereas 
u;^~^k) describes antiferromagnetic (AFM) one-magnon 
excitations. In the whole region occupied by the M = 1/2 
phase, the dispersion relations l(2Jl exhibit minima at the 
center (fc — 0) as well as at the boundary {k = tt) of the 
Brillouin zone. Close to these wave vectors the magnon 
energies take the generic form 
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Here A^ (h) and m^, are, respectively, the magnon 
gaps and effective masses. Since loq '{k) is a Goldstone 

mode, we have Aq (0) = 0. 

The magnon modes in Eq. (0J become gapless on dif- 
ferent parts of the phase boundary marked by the points 
a, b, c, and rf; uJq (k), ujir (fc), and ujir (k) are gap- 
less, respectively, on the pieces ab, be, and cd. At the 
special point b (analogous to the point E), both AFM 
modes Wg (k) and wJr (k) are gapless. Similar spe- 
cial points have been found in other ID spin models 
with folding excitation spectra.^ In the case when the 
related magnon masses in Eq. |0J diverge, such points 
produce cusps in the magnetization curves. Otherwise, 
as in our case, these are crossover points marking the 
boundary between regions with different low-lying exci- 
tations which belong to a single excitation branch. On 
the other hand, at the special point c with coordinates 
{J*,h*) = (0.75 - V2/6, 1.5 - \/2/6), two different types 
of modes - the FM mode wjr (fe) and the AFM mode 
i^TT (k) - reach the ground state. The point c belongs 
to the special line h = 0.75 -I- J2 on which the energies 
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of both magnon modes coincide: lu^ (k) — ujir (k)- Us- 
ing the picture of "particle-hole" excitations defined by 
^phik) = u}'^^\k) +tj'^^^(fc), we find that close enough 
to (J2,/i*) the dispersion relation ujph{k) takes the rela- 
tivistic form 



^jph{k) = y . 



Kn 



+ ^lH{k~^Y 



(5) 



While the velocity Vph remains finite at (J2 , h*), the ex- 
citation gap vanishes as lS.ph oc (J2 — J2) • Clearly, 



the discussed structure of the low-lying excitations in a 
vicinity of (J^j^*) remembers the known structure of 
low-lying excitations close to the commensurate Mott- 
insulator-superfluid transition,^ 

In fact, the outlined peculiarities of the (J2,/i) phase 
diagram - basically deduced from the linear SWT ap- 
proximation - share a number of common features with 
the {J/U,fi/U) phase diagram of the boson Hubbard 
model on ID periodic lattices. Before discussing this 
issue, let us present our quantitative estimates for the 
boundary positions based on higher-order spin-wave se- 
ries and numerical ED data. Unfortunately, already the 
first-order corrections in the SWT series exhibit diver- 
gences of the form 6^^^^, S being the distance in the 
parametric (J2,/i) space from the point (J^j/i*). Thus, 
reliable quantitative results can be achieved only for large 
enough S. As an example, we show in Fig. ^the first- 
order SWT result for the critical field he defined through 
the relation Aq ^ (he) = 0. In spite of the relatively large 
correction of the boundary position, we observe a good 
agreement with the ED data. In fact, one may expect fur- 
ther improvement of the theoretical result in a second- 
order SWT approximation!^ Relatively large finite-size 
effects in the ED data can be indicated only in the vicin- 
ity of the tip of the region occupied by the AI == f/2 
phase. It is remarkable that the peculiarities of the phase 
boundary - found in the framework of the linear SWT 
approximation - are reproduced by the ED data. In par- 
ticular, the cusp at point B - signaling a change in the 
type of low-lying excitations ~ can easily be indicated by 
the ED data. This is important because the ED data 
(due to finite-size effects) do not fix the position of the 
tip. 
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FIG. 2: Magnetization curves of the model Q for different 
values of the frustration parameter J2 {= 0.38, 0.42, 0.5, and 
1). 771 = M/Ms is the reduced magnetization, M being the 
magnetic moment per rung, and Ms = 3/2. The symbols 
mark the midpoints of the magnetization steps, as obtained 
from the ED data for L — 12. For clarity, we do not show the 
data points for J2 = 0.38 and 0.42 in the interval he < h < hs 
since they closely follow the presented data for J2 — 0.5. 



To compare with the phase diagram of the boson Hub- 
bard model, recall that the magnetic field h plays the 
role of a chemical potential and the magnetization Af is 
an analog of the conserved density of particles p. As is 
well known, the Luttinger liquid phase is an analog of 
the superfluid phase in the boson model, and the plateau 
phases can be viewed as Mott-insulator phases, -'^^ char- 
acterized by an integer p, a finite gap for particle-hole 
excitations, and zero compressibility, k = dp/dp=0. Ex- 
amples of magnetization curves m(h) for different param- 
eters J2 are presented in Fig. |21 These are analogs of the 
p{p) curves for the ID boson Hubbard model,-- Note 
that in the case J2 = 1 the curve m{h) describes the 
reduced magnetization of the spin-| antiferromagnetic 
Heisenberg chain»iS The magnetization curves also indi- 
cate a divergence of the zero- field magnetic susceptibility 
X = x(^2) as J2 approaches the isolated critical point Jc 
(J2 — > Jc ± 0), in accordance with the generic Luttinger 
liquid relation'^ x(/i, J2) = K{h, J2)/[TrVs{h, J2)] and the 
observation that the spin- wave velocity Us(0, J2) -^ 
as J2 ^ Jc ± Oi^ Here K{h, J2) is the Luttinger liquid 
parameter at finite h and J2. Up to now, only a few 
studies analyzing such ferromagnetic quantum critical 
points in ID have been published>i2*i^4i^ In particular, 
as recently pointed out by Sengupta and Kim^i^, in the 
isotropic case the ferromagnetic transition at Jc can be 
different from the ferromagnetic transition in the Ising 
limilii^ due to quartic terms in the effective action which 
couple the longitudinal mode to the transverse spin- wave 
modes. The magnetization curves also indicate that the 
canted magnetic phase is characterized by a finite sus- 
ceptibility x(0, J2)i as can be expected since this state 
may be viewed as a kind of ferromagnetic Luttinger liq- 
uid phase ii& 

Discussed instabilities at he and h^ share the same 
physics as those at the lower (hd) and upper (/ic2) crit- 
ical fields in Haldane-gap chains in a magnetic fieldii^ 
In particular, the instability at he may be viewed as 
a condensation of a dilute gas of magnon excitations, 
as in the case of Haldane-gap chains slightly above the 
lower critical field /ici >i24i2i2£ At the critical field he, the 

gap of the AFM mode uJq (k) vanishes and, as a re- 
sult, a macroscopic amount of AFM magnons condense 
onto the ground state. The collapse of quasiparticles is 
prevented by the repulsive on-site quasiparticle interac- 
tion. For h > he, the magnetization M continuously ex- 
ceeds the quantized value 1/2. The instabilities related 
to the FM modes, as in the case of a Haldane chain at the 
upper critical field hc2, are expected to share the same 
critical propertiesi2i In fact, the transitions discussed 
are of a commensurate-incommensurate type^SS like the 
insulator-superfluid transition in the ID boson Hubbard 
model on the phase boundaries where the density contin- 
uously changes from p = 1 to some incommensurate value 
p > 1 or p < 1. Since magnons behave as a dilute gas 
of hard-core bosons, in a first approximation the specific 
form of the interaction is irrelevant i22i2^ This results, in 
particular, in the universal behavior of the Haldane chain 



magnetizationiSiSi M = [2A//(/i - hci)]^^^ /(ttUs), where 
Ah is the Haldane gap and Vg (~ 2.49) is the spin-wave 
velocity. Similar universal expressions can be obtained 
close to the phase boundaries of the {J2,h) diagram of 
our model. Using the equivalence of ID hard-core bosons 
and free spinless fermions, this implies, in particular, the 
Luttinger liquid parameter iiT = 1 on the entire boundary 
of the M = 1/2 phase as well as on the /is( J2) boundary, 
excluding the special point at the tip of the region occu- 
pied by the M = 1/2 phase. Respectively, the spin-spin 
correlation exponents take the values rjz = 2 and ?7x = 2 ' 
and the magnetization exponent (3 = 1/2. Based on the 
established picture of the low-lying excitations close to 
the point (J|, /i*) and using the analogy with the boson 
phase diagram, it is plausible to assume that the latter 
critical point describes a Berezinskii-Kosterlitz-Thouless 
(BKT) type transition. Although our ED data do not fix 
the position of the tip, it is easy to indicate the typical 
pointed-type shape of the boundary^^ resulting from the 



specific BKT particle-hole excitation gap Ap/j close to 
this transition: Ap^ '^ exp {const j ^ J^ — J2). As might 
be expected, the SWT does not reproduce the correct 
behavior of the gap Ap/j close to this special point. 

In conclusion, we have demonstrated that a generic 
mixed-spin ladder exhibits a magnetic phase diagram 
which is closely related to the one of the boson Hub- 
bard model in one spatial dimension. Similar magnetic 
phase diagrams may be expected for arbitrary quantum 
spin numbers S\ and S2, provided that s\ 7^ S2 and 
5i + S2 =half-integer. 
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